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Based on a data sample of an integrated luminosity of 57.4 fb~ at the T(4S) resonance taken 
with the BABAR detector using the SLAC PEP-II asymmetric-energy B-Factory, hadronic decays 
of B meson with the final states D 1 -*' and three pions are studied. The study is performed by fully 
reconstructing the exclusive decays of B — > D'*^ai(1260) and the non-resonant modes B — > p°n 
and B — » £)'*' rnvn. The current status of the study is presented and a dominant B — > D^ai (1260) 
decay is shown. 

I. INTRODUCTION 

Hadronic B decays provide important information on both the weak and hadronic interactions of heavy flavored 
mesons. The quarks produced in such nonleptonic weak decays can arrange themselves into hadrons in many ways. 
The final state is linked to the initial state by QCD processes. The theoretical description of hadronic decays of heavy 
mesons invoke the factorization approximation and heavy quark effective theory (HQET) [1]. 

The dominant hadronic decay modes of the B meson involve tree-level diagrams where the b — > c transition leads to 
a charmed meson and an external W boson, which often emerges as a charged meson 7r, p or ai(1260). In the models 
based on the Bauer-Stech-Wirbel approach [1,2], two parameters a\ and a,2, which describe the QCD hard-gluon 
corrections for external and internal spectator processes respectively, provide important clues into the role played 
by the strong interaction in these two-body tree-level decays. Parameters a\ and ai are expected to be process 
dependent, but previous experimental data can be described with universal values (a\ « 1 and ai « 0.2). Recent 
results of B° -> D^°tt from CLEO and Belle (both indicate ai to be « 0.4) begin to show the expected process 
dependence of ai [3,4]. A precise measurement of B decaying to Z)(*'ai(1260) provides an interesting approach to 
estimate the constants a\ and ai and to test the factorization hypothesis. Since the QCD interaction between the 
product quarks continues after the weak transition takes place and after hadron formation, an understanding of final 
state interactions (FSI) is very important. The study of the non-resonant modes B — > p°ir and B — > D^nmr 
may provide important information on this mechanism. 

Hadronic B decays also play a significant role in CP violation study. A precise measurement of B° — » D^~af 
is important in measuring sin(2/5 + 7), where and 7 are the angles of the Unitary Triangle related to the CP 
asymmetries in B decay [5]. Presently B° — > D*~a^ is also suggested as a test of chirality [6]. 

Studies of B — » D^a\ and of their non-resonant decays were previously performed by CLEO and ARGUS [7,8,9,10]. 
Because of limited statistics, the uncertainties of the branching fraction measurements are relatively large [11], and 
some of the detailed substructure of the resonance is not available. BABAR has collected a large data sample of BB 
pairs at the PEP-II B- Factory [12], making possible an improved study of these modes. 

In this paper we present the status of our study. The BABAR detector and data are described in Section II, and the 
event selection with data analysis is presented in Section III. A summary is given in Section IV. Charge conjugation 
is always implied in this report. 

II. THE BABAR DETECTOR AND DATA SET 

The data used in this analysis consists of 62.2 million BB pairs, corresponding to an integrated luminosity of 
57.4 fb" 1 , collected with the BABAR detector at the T(45") resonance. The BABAR detector [13] was constructed 
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to observe CP violation and was motivated by the measurements of B mixing made almost two decades ago [14]. 
Charged particles are reconstructed with a five layer, double-sided silicon vertex tracker (SVT) [15] and a 40 layer drift 
chamber (DCH) with a helium-based gas mixture [16], placed in a 1.5 T solenoidal field produced by a superconducting 
magnet. The charged particle momentum resolution is approximately (Spr) 2 = (0.0013jjt) 2 + (0.0045) 2 , where pt is 
given in GeV/c. The SVT, with a typical single-hit resolution of 10 /um, provides measurement of impact parameters 
of charged particle tracks in both the plane transverse to the beam direction and along the beam. Charged particle 
types are identified from the ionization energy loss (dE/dx) measured in the DCH and SVT, and the Cherenkov 
device (DIRC) [17]. Photons are identified by a CsI(Tl) electromagnetic calorimeter (EMC) with an energy resolution 
a{E)/E = 0.023 • {E/GcY)- 1 ' 4 © 0.019 [18]. 

III. EVENT SELECTION AND DATA ANALYSIS 

B mesons are reconstructed in modes of B decaying to and three pions, which may form the resonance oi (1260) 
and non-resonant p°ir and 7T7r7r. They contain B° — ► D*~ af , B + — > D*°af , B° — ► D~af , B + — > D°af , and their 
corresponding non-resonant modes p ir and D^tttttt. 

To reconstruct the charmed mesons D and D* and the light mesons in the final states, the charged tracks are 
required to have a distance of closest approach within ±10 cm in z and 1.5 cm in radius of the average beam spot 
position, and at least 12 hits recorded in the DCH. Kaons and pions are selected, depending on their track's mo- 
mentum, based on the dE/dx information from the SVT and DCH, as well as the Cherenkov angle and the number 
of photons measured with the DIRC. For each detector component d = (SVT, DCH, DIRC), a likelihood L*f (LJ) 
is calculated given the kaon (pion) mass hypothesis. Photons are selected as showers in the calorimeter that have a 
lateral shower shape consistent with the expected pattern of energy deposits for an electromagnetic shower and with 
an energy E 7 > 30 McV. 

The D° meson is selected from the decays K~n + , K~n + ir° and K~TT + n~ir + , while the D + meson is selected from 
K~-k + tt + and K s n + modes. The track momentum is required to be greater than 200 MeV/c. ir°s are reconstructed 
from two photons with the sum of their energies E 11 greater than 200 MeV/c and an invariant mass 0.120 GeV/c 2 
< m 77 < 0.150 GeV/c 2 (about ±2.5 times the resolution of the 7r° invariant mass). A mass constraint fit is applied 
to 7T° candidates, which improves the energy resolution from 2.4% to 1.8%. Kg mesons are reconstructed from two 
pion tracks with opposite charge. The invariant mass of the candidate tt + it~ is required to be |m wvr — 497.7 MeV/c 2 
| < 15.0 MeV/c 2 and the momentum of the Kg greater than 200 MeV. In Kg selection the angle between the flight 
direction, which is the vector which points from the primary vertex to the 7r + 7r~ candidate's vertex, and the direction 
of the 7r + 7r~ candidate's three momentum is used as a cut. A vertex fit is applied, and a probability of x 2 greater than 
0.1% is required. The invariant mass of the candidate of D° — ► K ~n + , K ~7r + 7r~7r + , and D + — > K ~ir + n + , K s n + is 
required to be within 3cr from the true value of the mass. As the combinatorial background of the D — > _ftr~7r + 7r° is 
larger due to the presence of the n°, a 2a mass cut is applied on this mode. The tracks of the D meson are required to 
originate from the same point; therefore a vertex fit is applied with the requirement of a probability of \ 2 greater than 
0.1%. Then a combined vertex and mass constraint fit is applied and a convergent fit is required. The mass-constraint 
fit changes the momenta of the D meson daughters forcing the invariant mass to be the nominal one, thus improving 
the D energy and momentum resolutions. 

Reconstructed D°s are combined with soft pions tt~ (it ) to form D*~ (D*°) candidates. The soft pion momentum 
is required to be less than 450 MeV/c. The requirements of E 17 > 200 MeV for 7r° and the minimum momentum 
requirement of 200 MeV/c for the pion are removed for soft pions. D* candidates are selected by the requirement 
that the mass difference between D* and D°, Am = raw — m D a, lies within ±3cx (a is the resolution of Am) of 
the nominal mass difference. Then D* + — > D°ir + candidates are refitted with the beam-spot constraint to improve 
the angular resolution for the soft pion, and a convergent combined vertex and mass constraint fit is also applied. 
For D* Q — > D°tt°, no vertex fit is applied, but a kinematic and mass constraint fit are applied and required to be 
convergent. 

p° mesons are reconstructed from ir + ir~ pairs with the requirement of pion momenta greater than 200 MeV/c and 
the invariant mass satisfying \M n + n - — 0.770 GeV/c 2 | < 0.15 GeV/c 2 . af mesons are reconstructed by combining 
the selected p° and a pion with the pion momentum greater than 200 MeV/c and the invariant mass m p o n ± to be 
between 1.0 GeV/c 2 and 1.6 GeV/c 2 . A vertex fit is applied to the candidates of p° and af, and a \ 2 probability 
greater than 0.1% is required. 
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The B° and B + mesons are reconstructed by combining selected D or D* with a\ candidates, p°7r or tyutt. The 
variables AE and rriES are used to define the B signal. AE = E* B — Ej* eam is the energy difference between the 
energy of the B candidate and the beam energy in the T(AS) system, i.e., E* B is the center of mass energy of the 
B candidate and Et, 



is the center of mass beam energy. The B signal is expected to peak at \AE\ = 0. The 

2 where p* is the center 



beam energy substituted mass of B candidate, mssi is defined as mss — \J E'b eam 
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of mass momentum of the i-th daughter of the B candidate. The resolution in this variable is limited by the beam 
energy spread, which is about 2.7 MeV for BABAR. Since the final states involve many tracks and more than one 
candidate can be found in the event, the combinatorial background is high. To suppress such background, we chose 
for each mode the B candidates whose daughter D or D* masses are most consistent with their nominal given masses. 
Then among them the candidate with smallest \AE\ is selected. This selection will remove most of the combinatorial 
background and keep the best candidates. 



One major source of background comes from the qq continuum background. The continuum background is due to 
the large non-resonant fraction of the hadronic cross-section, approximately 75% at the T(4S I ), from direct e + e~ — > 
QQ (q — c 7 s i u j d). To select BB events from T(45*) data and reduce the continuum background, the ratio of second 
to zeroth Fox-Wolfram moment, R2, is used [19]. In the rest frame of T(4S*), R2 approaches zero for spherical events 
and one for jet-like events (Fig. 1). R2 < 0.4 is required in event selection since the momenta of B's are very small 
and BB events are spherical, whereas qq continuum events are jet-like in the frame of T(45). 




R2 

FIG. 1. Distribution of R-2- The solid-line dis- 
tribution is for BB Monte Carlo events and the 
dashed-line distribution is for continuum qq Monte 
Carlo events. 
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FIG. 2. Normalized distribution of |cos(#t)|. 
The hatched area is the distribution for BB Monte 
Carlo events. The solid-line and dashed-line his- 
tograms are for continuum cc and qq (q = u, d, s) 
Monte Carlo events respectively. (The vertical axis 
has arbitrary units). 



The thrust angle 8t, which is the angle between the thrust axis of the B candidate and the thrust axis of the remain- 
ing tracks in the event, is also used to suppress the continuum background. Monte Carlo simulation shows that the 
distribution of the thrust angle for continuum background events and of BB events are quite different (Fig. 2). The 
distribution of | cos(#t)| is flat for BB events and peaks at 1.0 for continuum qq events. We require | cos(9t)\ < 0.70 
for the B~ — ► D°a^[ , D°p°tt~ and D TT + Tr~n~ modes with \cos(9t)\ < 0.85 for all other modes. 



A. Resonant mode B w -> D^ai (1260) 

In this decay B mesons are selected in the modes B° — > D*~a^, B + — ► D* a a^, B° — > D~al and B + — > _D°a^. 
The meson is a very broad (r k, 400 MeV) isovector (1 = 1) state with / = 1 orbital excitation and J p = 1 + . 
Fig. 3 shows the Monte Carlo simulated spectra of mass and momenta in the decay of B — > Z)(*)ai(1260). arc 
reconstructed from a combination of selected p° and a charged pion, with the pion momentum greater than 200 MeV/c 
and the invariant mass m p o^± between 1.0 GeV/c 2 and 1.6 GcV/c 2 consistent with the mass. An additional cut on 
the center of mass momentum of the a\ candidate is applied with p* ai > 0.5 GeV/ c. A vertex fit is performed and a x 2 
probability greater than 0.1% is required. All selected D° and D + candidates from B decays or D* decays are required 
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to have a momentum p* D in the T(45) frame greater than 1.3 GeV/c. The AE distributions for B — > Z)(*)ai(1260) 
with m£s > 5.27 GeV/c 2 are shown in Fig. 4, while the rngg distributions with \AE\ < 2.5ctae are shown in Fig. 
5 and Fig. 6. 
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FIG. 3. Monte Carlo simulated spectra (normalized distribution with arbitrary unit invertical axis) for (a) a{ mass and (b) 
momenta in the decay of B — » D^ai (1260), and (c) the mass of p° from oi. 



The ities background can be separated into continuum and BB components. The BB background component is 
the result of mis-reconstructing other BB decays. The relative contributions and the overall amount of background 
varies decay-mode by decay-mode depending primarily on the multiplicity of the B decay. We use generic Monte Carlo 
BB data with the signal modes B — ► D^a\ and non-resonant modes B — > D^p°n and B — ► D^tttttt removed, and 
continuum Monte Carlo qq data to model the background and find that the background shape of data can be well 
characterized by Monte Carlo (Fig. 7). The Monte Carlo backgrounds are fitted to an Argus function [20] and the 
obtained shape parameters of such functions are used in the fitting of tues of data as shown in Fig. 5 and Fig. 6. 
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FIG. 4. AE distributions for B° -> D ( *^a+ with m E s > 5.27 MeV/c 2 . 

In the signal region of B — > D^a\, some events of the non-resonant decays of B — ► D^p°tt and B — > D^tvtttt 
can pass the event selection criteria of B — > D^a\ and form a peaking background. This is confirmed by BB Monte 
Carlo data. Fig. 8 shows that peaking background of all the non-resonant modes underlays the Monte Carlo signal 
of B — > D^a\. The estimation of contributions from non-resonant modes to B — > D^ai modes will be given in 
the next subsection. It is also possible that the decay of B° — > D*^tt + tt + may pollute the signal of B~ — ► Z?°a^ by 
misreconstructing p° from one direct tt + and the slow tt~ from Z)* - , where the mass of p a ir + lies in the a\ region. 
Monte Carlo study shows that compared with the decays B — > D*p°-K~ and B —> D*ir~Tr + ir~, its contribution is 
quite small and can be neglected. Another possible source of background is from the decay B~ — > D°(2420)7r~ (p~) 
and B~ L>S?(2460)tt- (p"), where or D% decays to D*+tt- [21]. Monte Carlo study also shows that their 
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contributions are quite small and can be neglected. 
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FIG. 5. m BS distributions for B° -> £> ( * )_ a+ with |AE| < 2.5<tab. 
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FIG. 6. m ES distributions for B~ -> D ( * )0 aj" with \AE\ < 2.5a AE . 
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FIG. 7. rriES distributions of B — * D^ai from data (dots) and background modeled by Monte Carlo (shaded area). 

The event selection efficiencies for each submodc have been studied using signal Monte Carlo data and are reported 
in Table 1. Estimates of the statistical and systematic error of branching fraction measurements are also listed in 
the table. The systematic error is mainly considered from the following sources (1) the branching fraction errors of 
intermediate decays, (2) the error of the integrated luminosity (about 1.1%), (3) the error on the acceptance efficiency 
due to Monte Carlo statistics, and (4) the error on the charged track reconstruction efficiency, which is 1.2% per 
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charged track. An additional error of 1.6% is added in quadrature to account for the uncertainty in the soft pion 
reconstruction efficiency. 



TABLE I. Event selection efficiencies and branching fraction error estimation. 
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FIG. 8. Peaking background of non-resonant modes (dotted/colored line) underlays the signal of B — > D^ai (BB Monte 
Carlo data, with m„o v < 4.0 GeV/c 2 ). 



B. Non-resonant modes B — ♦ D'*'p°7r and B — > D^-k-k-k 

For the decays of B — > p°ir and £? — > D 1 -*) it-kit, the invariant masses of p°7r and 7T7r7r are required to be less 
than 4.0 GeV/c 2 and center of mass momenta of p°ir and tttttt greater than 0.5 GeV/c. The tubs distributions for 
B — ► D^p a Tr are shown in Fig. 9. 
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The B signal region of B -» D^pPir, which is |to£;s - 5.280 GeV/c 2 | < 3cr mfiS and |A£| < 2.5a A E (see Fig. 
9) also contains B — > D^-*'ai and £> — > D^tttttt candidates. Monte Carlo study shows that the kitty invariant mass 
spectra for B — > £>Wp°7r~ and 5 — > dW-tt^tt+tt - are different. But when the cuts \m„+ v - - 0.77 GeV/c 2 | < 0.15 
GeV/c 2 and center of mass momentum of ir~ir + Tr~ greater than 0.5 GeV/c are applied to both modes, the ir~ir + 7r~ 
mass spectra for B — > £K*)p°7r - and i? — > £>(*'tt~7t + 7t~ are almost the same, as shown in Fig. 10. So, we can use the 
7r~7r + 7r~ spectra of B — ► D^*'p°n to represent the spectra of I? — > n~ tt + n~ . Monte Carlo also shows that the 
7r~7r+7r~ spectra of modes B° — > D~p Q n + , B° — ► D*~p Q n + , B~ — > D°p n~ and B~ — > D* p n~ are quite similar 
because they have similar topology and kinematics (see Fig. 11). 
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FIG. 10. Monte Carlo simulated distributions of -k~tt + tt~ mass for B -> D^p°TT~ (upper left), B -> £> w 7T 7r + 7r (upper 
right) and 7r _ 7r + 7r _ mass (lower left) for B -> D (t) p°-K~ (solid-line) and B -> 7T 7r + 7r (dotted/colored-line) with the cuts 
Im^+jr- —0.77 GeV/c 2 | < 0.15 GeV/c 2 and center of mass momentum p* K - 1T+7r - > 0.5 GeV/c. (Normalized distributions with 
arbitrary unit in the vertical axis). 
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FIG. 11. Monte Carlo simulated p°n + mass distribution for decay B° -> D~p°n + , B° -» D*~ p°ir + , B~ -> D° p°n~ , and 
B~ — > D*°p°tt~ . (Normalized distributions with arbitrary unit in the vertical axis). 



The invariant mass distributions of p°7r + in the i? signal region are shown in Fig. 12. The shaded areas in the 
figures are the distributions from vtes sideband region (5.22 GeV/c 2 < ities < 5.26 GeV/c 2 ), and were scaled with 
the ratio of the integral of the Argus function in the rriES sideband region to the integral in the ities signal region. 



To estimate the fraction of ai(1260) in the p°ir~ mass distribution, the p°7r _ mass distribution, with the ities 
sideband background subtracted, are fitted to a Breit-Wigner distribution as ai(1260) signal plus the distribution of 
non-resonant p°ir + mass. The fits of B° modes are shown in Fig. 13. The non-resonant p°ir + mass distribution are 
from the reconstruction of the p°Tr + mass in signal B — > £>(*'p°7r + Monte Carlo data. It is found that the Oi(1260) 
signal is dominant in the mass spectra of p°ir + in both B° and B~ decay; therefore, the decay of B — > £>(*)ai(1260) is 
dominant in B decay with final state D^*> p n. To separate the non-resonant mode B — > D^p°n and B — » D^tttttt, 
a Dalitz plot analysis will be important. 
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FIG. 12. p°7r mass distribution in the tues signal region with the shaded area showing the ities sideband background. 
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IV. SUMMARY 

The hadronic decay of B — > Z)(*)ai(1260) and non-resonant modes B — > D^p°n, B — > D^nnn are studied with a 
large data sample. The preliminary result shows that i? — > £>^*- ) ai(1260) is dominant in _B decays with the final state 
D^p°tt. The errors on the branching fraction measurements are estimated and improved compared with previous 
measurements. Further study is ongoing and final results will be published soon. 
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